Molecular adsorption in graphene with divacancy defects 
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We have investigated theoretically the adsorption of molecules onto graphene with divacancy 
defects. Using ab-initio density functional calculations, we have found that O2, CO, N2, B2 and 
H2O molecules all interact strongly with a divacancy in a graphene layer. Along with a complex 
geometry of the molecule-graphene bonding, metallic behavior of the graphene layer in presence of 
CO and N2 molecules have been found with a large density of states in the vicinity of the Fermi level 
suggesting an increase in the conductivity. The adsorption of N2 is particularly interesting since the 
N atoms dissociate in the vicinity of the defects, and take the place where the missing C atoms of 
the divacancy used to sit. In this way, the defected graphene structure is healed geometrically, and 
at the same time doped with electron states. 

PACS numbers: 73.20.Hb, 81.05.Uw, 71.15.Nc 



I. INTRODUCTION 

In the last few years, graphene has been considered 
as one the most sensational materials for its many ex- 
otic properties.^ Graphene is known theoretically for sev- 
eral decades but the real breakthrough came when it 
was discovered experimentally in 2004.^2! The electronic 
properties of this sp^-bonded two-dimensional allotrope 
of carbon have striking similarities with those of mass- 
less relativistic Dirac fermions. The linear band disper- 
sion at the so called Dirac point is a special feature of 
the graphene band structure. Apart from this, graphene 
has been considered as an excellent candidate for future 
electronics due to the extraordinarily high mobility of its 
charge carriers. As today's Si based electronics faces the 
fundamental limitation of spatial dimension lower than 
a few nanometers, graphene may stand as a potential 
alternative. 

The electronic structure of graphene may be engi- 
neered by the interaction with supporting substrates, 
gate voltag^ or by chemical doping with molecules. Re- 
search on graphene is being pursued in different direc- 
tions to achieve novel properties. One of the routes is to 
tune the properties of graphene by controlled defect for- 
mation through chemical reactions. In a recent papeii^, 
it has been shown that chemically treated graphene 
nanosheets have oxygen-containing defects, which were 
identified by X-ray absorption spectroscopy (XAS). The 
oxygen containing groups were introduced by the acid 
treatment of the nanosheets, and follow a reaction path- 
way similar to what happens for the case of carbon nan- 
otubes with the addition of chemical reagents to car- 
bon network defects in the graphene sheets. Electronic 
structure calculations revealed that the defect states arise 
from a network of vacancy defects, where carbon atoms 



have been removed. Moreover, it was shown theoretically 
that the carbon atoms around such a lattice defect show 
metallic behavior due to the presence of high density of 
states at Fermi level in contrast to a vanishing density of 
states (DOS) for a pure graphene sheet. Such a defect 
would also be very reactive due to the highly strained 
structure, and hence efficiently undergo acid-catalyzed 
water addition, the initial product of which rapidly rear- 
ranges to form a pyrane-like structure that possess the 
characteristics found in the XAS of acid-treated carbon 
nanosheets (CNS). 

One effective route to induce charge carriers in 
graphene would be by molecular doping, e.g., by NH3, 
NO2 etc. It has been shown by Hall measurements that 
NO2 introduces holes whereas electrons are introduced 
by doping with NH3. This way has opened up an ex- 
cellent playground for both experimentalists as well as 
theoreticians. Wehling et al.'^have demonstrated by den- 
sity functional calculations that a NO2 molecule acts as a 
strong acceptor and produces special features in the den- 
sity of states (DOS) that is suitable for chemical sensor 
applications. It was also shown recently^ by density func- 
tional calculations that chemical functionalization of a bi- 
layer graphene can open up a gap up to 3 eV. It is note- 
worthy to mention that the detection of gas molecules 
adsorbed in graphene has been possible recently in ex- 
periments by Schedin et al^. These graphene based gas 
sensors have an extraordinary sensitivity to detect the 
change in local charge carrier concentration by minute 
amount due to the extremely low noise in graphene. Be- 
sides the characterization and modification of the prop- 
erties of infinite graphene sheets, a considerable effort 
has been made recently to predict properties of graphene 
nanoribbons (GNRs). Modification of electronic proper- 
ties using edge functionalization, substitutional doping 
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and also the transport and magnetic properties of GNRs 
have been studied in great detail.'^ The edges of GNRs 
resemble the vacancy defects in graphene and hence are 
relevant for the present discussion. 

The experimental information on the acid-catalyzed 
addition of water to edges and strained defects without 
dangling bonds, such as the thermally mobile 5/7 mem- 
bered rings are available from the research on carbon 
nanotubes (CNTs) and also some transmission electron 
microscope (TEM) work on graphites. The addition of 
water is reversible and heating at low pressure removes 
quite a lot of water molecules to form a very reactive sur- 
face, e.g., defects with dangling bonds that would relax to 
still quite reactive but more stable 5-8-5 or related lattice 
defects. The defects in carbon-bond-network has been 
studied in detail theoretically by Dias.^ The purpose of 
this paper is to study the chemical interaction between 
divacancy defect centers in graphene and chemisorbed 
molecules with the aid of theoretical analysis. Using den- 
sity functional based electronic structure calculations, we 
have analyzed the reaction between water molecules and 
defect centers and consequently the passivation of these 
highly reactive centers by several molecular species. 



II. COMPUTATIONAL DETAILS 

Here, we have performed ab-initio electronic structure 
calculations based on the density functional theory.'J^ 
The Vienna Ab-initio Simulation Package (VASPjfSl em- 
ploying a plane wave basis was used within the all elec- 
tron projector augmented wave (PAW) method. The 
generalized gradient approximation (GGA) was used as 
the exchange-correlation functional. GGA reproduces 
very well the C-C bond length (1.42 A) in graphene. 
The kinetic energy cut-off was chosen to be 400 eV. We 
have optimized the geometries by reducing the Hellman- 
Feynman forces down to 0.01 eV/ A. 

We have considered five different molecules and studied 
their interactions with a divacancy center (e.g. as shown 
by the C network in Fig. 2a) in a graphene sheet. The 
molecular species we have considered are O2, N2, B2, CO 
and H2O molecules. The graphene sheet was modeled 
with a 8x8 supercell in the lateral direction with a 15 
A vacuum in the vertical direction. This is a relatively 
large lateral supercell to reduce the interaction of the 
adsorbed molecules with their periodic images and can 
be approximated well by isolated ones. 

From the total energy calculations, we have extracted 
the chemisorption energies for all the molecules consid- 
ered. The chemisorption energy of the molecule is calcu- 
lated as 
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where Euv+m, Euv and Em denote the total energies 
of a graphene sheet with a divacancy and the molecule, 
a graphene sheet with a divacancy, and a free molecule, 
respectively. The calculated values of E'^'^ are -8.44 eV, 



-4.53 eV, -13.83 eV and -3.86 eV respectively for O2, 
N2, B2 and CO molecules. All these numbers are quite 
large and they indicate a very strong binding between 
the molecules and graphene divacancies. A large bind- 
ing energy between O2 molecule and defected graphene 
indicates it is rather easy to oxidize defected graphene. 
In case of B2 , a very large value may arise from the fact 
that this molecule doesn't form naturally as a diatomic 
molecule. 



III. RESULTS AND DISCUSSIONS 



A. Geometry 



(a) 




FIG. 1: (Color online) (a) Initial structure of II2O, (b) final 
structure after geometry optimization. C, O, and H atoms are 
shown in yellow (light gray), red (dark gray) and turquoise 
(medium gray) colors respectively. 



Analysis by chemical models suggests that a defect 
lacking two carbon atoms forming a 5-8-5 lattice defect is 
more probable than a single vacancy. These defect cen- 
ters are highly reactive as the carbon atoms at the de- 
fect edges have dangling bonds. When this non-perfect 
induced graphene is treated with water, chemical rear- 
rangement occurs as shown in Fig. [T| In our study we 
have used a II2O molecule as shown in Fig. [ija) in the 
beginning of the simulation. In the presence of H2O, 
one C atom at one side of the defect is saturated by H 
and the other C atom is saturated by OH. The two C 
atoms at the remaining pentagon of the lattice defect 
are left untouched. This is one of the probable configu- 
rations existing in the reaction. The final configuration 
after geometry optimization is shown in Fig. [ijb). The 
initial quite strained phenolic product is rearranged to 
a more relaxed pyrane-type structure, with the O atom 
in a position that produces a hexagonal pattern and the 
two hydrogen atoms saturating the remaining two carbon 
atoms of the original defect. The H atoms are aligned at 
two opposite faces of the graphene sheet making an angle 
of around 28° to the graphene plane. The final structure 
is very stable, as all the dangling bonds are saturated. 

Similar to the case of H2O, a graphene with unsatu- 
rated dangling bonds will immediately absorb any molec- 
ular center close to it. We have considered O2, CO, B2 



3 



and N2 molecules interacting with the lattice defect in 
the C network of graphene. Fig. [2|a) shows the initial 
configuration of the simulation. The molecule is kept at 
a horizontal position above the vacancy region and the 
geometry of the system was optimized by minimizing the 
forces. The final optimized configurations are shown in 
Fig. [gb), [gc), gd) and IJe) for O2, CO, N2 and B2 
molecules respectively. 




FIG. 2: (Color online) (a) Initial structure for O2, CO, N2 
and B2 molecules. Final relaxed structures for (b) O2, (c) 
CO, (d) N2 and (e) B2 molecules. In all panels, C atoms are 
shown in yellow (light shade). In (b), O atoms are shown in 
red (dark gray), in (c), C and O atoms of CO are in gray 
(gray) and red (dark gray) respectively, in (d), N atoms are 
in turquoise (medium gray) and in (e), B atoms are shown in 
blue (blackish gray). 



The three molecular species behave differently while 
being adsorbed in the graphene sheet. The O2 molecule 
gets dissociated into two O atoms bound to two pairs of 
C atoms around the divacancy. The O atoms are moved 
in opposite directions out of the graphene plane, a situ- 
ation similar to that of H atoms described in Fig. [T] For 
the case of CO, the C and O atoms still remain bonded 
to each other, while being bonded to two pairs of carbon 
atoms. Moreover, the CO molecule doesn't remain in the 
plane of graphene. The C atom of the CO molecule stays 



in the plane making a usual hexagonal planar structure 
while the O atom moves a little bit out of the plane. 
The N2 molecule is completely adsorbed in the plane 
of graphene. The two N atoms are substituted in the 
normal positions of a graphene network. So, a graphene 
sheet with substitutional N atoms is achieved in this pro- 
cess. This may serve as a route to dope N atoms in the 
graphene network and hence introduce extra charge car- 
riers to achieve novel transport properties. In case of B2, 
the B atoms are buckled out of the graphene plane by 
about 0.55 A. 



B. Electronic structure 




E - Ep (eV) 

FIG. 3: (Color online) DOS for CO on graphene. In the 
panels. Cm, Om, Co (O) and Cg (C) indicate DOSs for C 
atom in CO, O atom in CO, C atom in the graphene sheet 
connected to O and C atom in the graphene sheet connected 
to C atom of CO molecule respectively. For all the cases, 
DOSs projected on different p orbitals are shown. 



We now present the electronic structure of graphene 
with CO, N2 and B2 chemisorbed to the lattice defect.. 
The case for the CO molecule is shown in Fig. |3] Gmi 
Om, Cg (O) and Cg (C) indicate the C atom in CO, 
the O atom in CO, the C atom in the graphene sheet 
connected to O and the C atom in the graphene sheet 
connected to the C atom of CO molecule, respectively. 
Projected DOSs on different atoms and orbitals show 
different behaviors. The Pz orbital-projected DOSs of C 
(Cm) and O (Om) atoms in the CO molecule and the C 
(Cg(0)) atom in the graphene sheet connected to the O 
atom of CO, have a rather sharp peak at the Fermi level 
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(Ei?). This feature is very similar to the mid-gap states 
observeci^ in graphene with defects. The C (Cg (C)) 
atom in the graphene layer connected to the C atom of 
CO has a usual graphene-like behavior, i.e., a vanishing 
DOS at Fip. For this atom the electronic structure is ac- 
tually quite similar to that of defect free graphene, with 
the Px and py orbitals being almost degenerate form- 
ing an occupied a state, whereas the Pz orbital forms a 
TT and TT* complex just below and just above the Fermi 
level, respectively. This picture is somewhat different for 
the other DOS curves in Fig.l, and it is clear that for 
Cm and Cg(0) the p^; and Pj^ orbitals are no longer de- 
generate, since symmetry is lowered due to the geometric 
reconstructions. 
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FIG. 4: (Color online) DOSs for (top left) H2O, (top right) 
O2, (bottom left) N2 and (bottom right) B2 molecules are 
shown. In the top left panel, C-p(O) and C-p(H) indicate 
the p orbitals of the C atoms of graphene bonded to O and 
H respectively. For O2, N2 and B2, C-p indicates the p or- 
bital projected DOS for one of the C atoms connected to the 
molecule. 

In Fig. m we have shown the atom and orbital pro- 
jected DOSs for N2, B2, H2O and O2 molecular species 
adsorbed on a divacancy in graphene. For N2 adsorption, 
we observe a large DOS very close to E^. In this case, E^ 
lies almost at a very sharp peak primarily of N pz charac- 
ter. One might suspect that this could lead to a Stoner 
instability and a spin-polarized ground state. In order 
to investigate this we also performed spin-polarized cal- 
culations. However, in these calculations the magnetic 
moment vanishes to result in a spin-degenerate state. 
The results of N2 adsorbed are nevertheless very inter- 
esting since the N molecules, after geometry optimiza- 
tion are found to take positions to generate a perfect 
honeycomb lattice. Hence, after N2 adsorption the diva- 
cancy defect has healed, generating a perfect graphene 



geometry, albeit with N atoms replacing some of the C 
atoms. This naturally provides electron doping into the 
lattice, and our results may hence provide a recipe for 
doping graphene with electrons. One may envision cre- 
ating divacancy (or similar) defects of graphene by ion- 
or electron-irradiation under vacuum conditions, and af- 
terwards flush the sample with N2 gas, and in this way 
provide a graphene like material, with electrons doped 
into the conduction band. In case of B2, a very sharp res- 
onance of Px character is observed just below the Fermi 
energy followed by a smaller peak of p^ character. The 
large DOS in the vicinity of Fermi level is absent for 
the cases of H2O and O2. In these pseudo- 
gap at the Fermi level is produced similar to the case of 
pure graphene. The DOS of the H2O adsorbed system is 
particularly interesting having in mind the experimental 
study of Ref , 4 , In the present study, peaks in the C and 
O atom projected DOSs are observed around 4 eV above 
the Fermi level. In the paper mentioned above, a peak in 
the XAS spectra was observed around 3-4 eV above the 
threshold. The occurrence of this peak was argued to be 
due to the presence of vacancy defects. Here, a similar 
observation is made in presence of the pyrane-type satu- 
rated lattice defect. One may conclude that unsaturated 
dangling bond vacancy defects as well as saturated bond 
vacancies saturated by reaction with water may be the 
cause for the modification of the electronic structure of 
graphene sheets. 



IV. CONCLUSIONS 

In conclusion, we have studied here the adsorption 
of different molecular species to a graphene with a di- 
vacancy, correpsponding to a highly likely 5-8-5 lattice 
defect. The calculated chemisorption energies for the 
molecules are quite large, with a magnitude which is of 
order 3-13 eV. Our ab-initio density functional calcula- 
tions indicate that N2 is adsorbed in the plane of the 
graphene network as substitutional impurities, whereas 
H2O, CO and O2 have complex geometrical constructions 
at the divacancy site. The B2 molecule more or less re- 
tains its bond length as in the isolated molecular state 
and is strongly bonded to the divacancy center in a po- 
sition a little bit off the graphene plane. The significant 
observation is that the inclusion of N2 in the graphene 
plane produces a perfect honeycomb lattice with elec- 
trons doped into the conduction band. A recipe for doing 
this experimentally is suggested. 
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